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Abstract 

This paper is concerned with experiments which measure CMB anisotropics 
on small angular scales. A certain coverage, a beam structure and a level of 
uncorrelated noise define each experiment. We focus our atention on the rever- 
sion of the beam average. In each experiment, we look for the best pixeliza- 
tion for reversion, namely, for the pixelization that -after reversion- leads to 
good maps containing right spectra for the most wide range of angular scales. 
Squared pixels having different sizes "smaller" than the beam radius iP FWHM ) 
are considered. For a given size, the following question arises: How well can we 
assign a temperature to each pixel? Various mathematical methods are used to 
show that, in practice, this assignation -beam reversion or deconvolution- only 
leads to right spectra for pixel sizes greater than a certain lower limit close to 



9 FWHM /2. This limit is estimated for negligible and relevant levels of noise and 
also for spherically symmetric and asymmetric beams. After this general study, 
we focus our attention on two feasible detectors (which have been proposed to 
be on board of PLANCK satellite). For each of them, we estimate the size of 
the most appropriate pixelization compatible with beam reversion, difraction, 
observational strategy et cetera and, then, we answer the following question: 
Which is the part of the angular power spectrum which can be extracted from 
appropriately pixelized maps after deconvolution? 

Subject headings: cosmic microwave background — cosmology:theory — large-scale 
structure of the universe — methods:numerical 
PACS: 98.70.Vc, 95.75.-z, 95.75.pq, 95.75. St 
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1 INTRODUCTION 

The main goal of this paper is to present a detailed analysis of the following problem: 
Given a beam structure, a level of noise, a certain partial coverage, and a pixelization, 
how well can we assign a temperature to each pixel?. In other words, how well 
can we deconvolve the beam to get appropriate temperatures in the pixels?. Here, 
temperatures are considered to be appropriate when the resulting maps lead to good 
physical spectra. Hereafter, any beam reversion leading to right spectra is referred 
to as a " S-deconvolution" . In the absence of noise, the possibility of performing a 
good S-deconvolution essentially depends on the ratio between the beam area and 
that of the chosen pixel. In practice, S-deconvolution is not feasible for too large 
values of this ratio; in other words, if we fix the beam, S-deconvolution is not feasible 
for too small pixel sizes. For a given beam and a certain mathematical method, 
there is a minimum pixel size allowing S-deconvolution. For values smaller than this 
minimum, too many pixels can be placed inside the beam and S-deconvolution is 
not possible. The minimum size corresponding to two S-deconvolution methods has 
been estimated in various cases. Both methods lead to similar minimum values of 
the pixel size around 9 FWHM /2. These values depend on the level of the uncorrelated 
noise produced by the instruments. 

It is worthwhile to emphasize that we are not interested in assigning temperatures 
to hundreds of pixels located inside the beam. This assignation can be useful in other 
contexts; however, in our case, the important point is that the spectra contained in 
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the S-deconvolved maps must be similar to the true physical spectra (up to the scales 
corresponding to the pixel size). Unfortunately, this type of deconvolution requires a 
moderate number of pixels inside the beam. This number will be estimated below in 
various cases. 

From § 2 to § 4, pixelization is considered in the framework of pure beam rever- 
sion, without analyzing particular experiments; however, in § 5, we focus our attention 
on PLANCK multifrequency observations and, then, pixelization is discussed taking 
into account both previous conclusions about beam reversion and some physical con- 
straints due to difraction, observational strategy, et cetera. 

In Saez, Holtmann & Smoot (1996) and Saez & Arnau (1997), the modified power 
spectrum 

327T 3 ramax 

E e( a ) = T7r B , 1U / C a (a)Pe(cosa) sinada (1) 

[M + X) Ja min 

was described. Functions Pg are the Legendre polinomial normalized as follows: 
/ PgPi'd^cosO) = [(2£+l)/87r 2 ]<5«/. As explained in those papers, this type of spectrum 
can be easily found from both theory and maps. Comparisons of the modified spectra 
obtained from theory with those extracted from simulated or observational maps are 
appropriate to take into account pixelization, partial coverage and beam features, 
simultaneously. Sometimes, the estimation and use of the well known Cg coefficients 
-although possible- is not the best procedure. In Eq. (P, the effect of pixelization is 
simulated by the angle a min , which is the angle separating two neighbouring nodes, 
while the angle a max depends on the area of the covered region; this angle is to be 
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experimentally obtained (Saez & Arnau, 1997), it is smaller than the size of the map 
and large enough to include as much scales as possible. The autocorrelation function 
C c {a) is 

C ^ = ({y) (^)) > ( 2 ) 

where a is the angle formed by the unit vectors Hi and n2, the angular brackets stand 
for an average over many full realizations of the CMB sky and, quantity (ST/T) a (n) 
is the temperature contrast in the direction n after smoothing with a Gaussian beam 
having a certain a = 0A259 FWHM . The modified spectra are used below to analyze 
some simulated maps. 

2 SIMULATIONS 

The angular power spectrum Cg = Em=-< \ a em\ 2 ' I (2£+ 1) is only an auxiliary element 
in our estimations. We are not particularly interested in any choice and, consequently, 
we have used the same spectrum as in Saez, Holtmann & Smoot (1996). It corresponds 
to the minimum cold dark matter model with a baryonic density parameter Q B = 0.03 
and a reduced Hubble constant h = 0.5. The Cg coefficients have been taken from 
Sugiyama (1995) and renormalized according to the four-year COBE data {Q rm s P s — 
1S/J.K , Gorski et al. 1996). Our simulations are performed by using the Fast Fourier 
Transform (see Saez, Holtmann & Smoot 1996 and Bond & Efstathiou 1987) and, 
then, a certain beam is used to average temperatures; thus, we obtain maps which 
must be deconvolved with the same beam. After S-deconvolution, the resulting map 
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must be compared with the initial one. 



3 IDEAL BEAM S-DECONVOLUTION 

Two methods are proposed to perform a S-deconvolution of the beam in the absence 
of noise (ideal case). The efficiency of these methods is verified and the limits for 
their application are discussed. For appropriate coverages and beam structures, the 
size of the smallest pixels compatible with beam reversion is estimated in each case. 
The conclusions obtained in this section are important to understand realistic S- 
deconvolution in noisy maps (§ 4). 

3.1 BEAM 

We begin with a Gaussian spherically symmetric beam. If the direction of the beam 
center is n, the measured temperature T(n) is given by the following average: 

T(n) = -^r / T*(9, <p)e- [e ' {R)]2/2a2 sin 9d9d<j> , (3) 
2ti<j z J 

where a defines the beam size, the angles 9 and <fi are the spherical coordinates of 
a certain pixel, the element of solid angle is dVL = sm9d9dcf) and, quantity 9* is the 
angle formed by the direction (9, 0) and the observation direction n. 

Small pixels can be considered as surface elements and, consequently, Eq. (|) can 
be discretized as follows: 

T(n) = J_YT*e-W )2/2a2 ^- U) 
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where the subscript i stands for the i-ih pixel; here, T* and dSi are the temperature 
and the area of the i-th pixel, respectively, and 9* is the angle formed by this pixel 
and the beam center. The exponential tends rapidly to zero as 9* increases beyond 
a; hence, only a small number of pixels are significant in order to estimate T(n). 
Furthermore, due to technical reasons, the beam could receive energy from a reduced 
number of pixels (not from all the significant pixels in an ideal infinite Gaussian 
beam). By these reasons, we assume that only q x q pixels are relevant and we give 
various values to number q; these pixels cover a square patch centered at the same 
point as the beam. 

An asymmetric beam of the form 

w = 1 c \aHe-e'?+a-H<t>-4>') 2 }/2^ (5) 
2na 2 K ' 

has been also considered. The parameter a defines the degree of asymmetry. 

Figure 1 illustrates, for q = 7, three situations corresponding to beams and pix- 
elizations considered below. 

3.2 COVERAGE 

Our choice of an appropriate partial coverage is based on some results obtained in 
previous papers. In Saez, Holtmann & Smoot (1996), it was shown that CMB maps 
close to 20° x 20° can be simulated -with good accuracy- neglecting curvature and 
using the Fourier transform. The effects of partial coverage were analyzed in detail 
in Saez & Arnau (1997), these authors proved that, although a 20° x 20° map does 
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not suffice to get the angular power spectrum, a few tens of 20° x 20° maps allow 
us to find the power spectrum with a high accuracy for t > 200. On account of this 
previous result, forty 20° x 20° maps are used in this paper to get each angular power 
spectrum. These maps would cover about 40% of the sky. In the above papers, it 
was also argued that uncorrelated noise is not expected to be important for coverages 
greater than or equal to 20° x 20° and for the noise level expected in modern CMB 
experiments. This is true in the sense that the power spectrum can be calculated in 
the presence of this noise, but the noise can be problematic for S-deconvolution as 
we will discuss below. These comments point out the interest of considering 20° x 20° 
maps and motivate our choice of these regions to begin with our analysis of the 
S-deconvolution procedure. 

3.3 PIXELIZATION AND EQUATIONS TO BE SOLVED 

In order to compute the integral in (|3|) using Eq. (U), the area dSi is not required 
to be independent on i; namely, no equal area distributions of pixels are necessary. 
Furthermore, if we take a large enough number of small pixels covering all the region 
contributing significantly to the integral (|3|), moderated variations in the pixel shapes 
are also admissible. In spite of these comments, equal area and equal shape pixeliza- 
tions are advantageous -at least from the mathematical point of view- as it is shown 
along the paper. 

Since a 20° x 20° region is approximately flat, small squares with edges of angular 
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lenght A6 = A0 = A define an approximately regular pixelization. The number of 
pixels per edge is N = 20/A, where angle A must be given in degrees. Our mea- 
surements would cover the pixelized region with a certain strategy. The beam center 
should point towards each pixel a various times and, then, the mean of the result- 
ing measures could be considered as the smoothed temperature T a corresponding to 
pixel a. From these T a values, a new temperature which does not involve the beam 
effect, T*, must be assigned to each pixel %\ namely, the beam must be deconvolved. 
All along § 3, it is assumed that sistematic errors have been corrected and that the 
uncorrelated instrument noise is negligible. If T a is the temperature measured by 
the instrument when the beam center is pointing towards the center of the pixel a, 
according to Eq. (^), we can write 

T a = £ A ai T* (6) 

i 

with 

A ai = J-e-WK*-)-)'/*^ , (7) 
2na 2 R 2 ' y J 

where T* is the true temperature in the i-th pixel; namely, the S-deconvolved tem- 
perature we are looking for. Similar equations hold for the asymmteric beam (El). 

It is evident that the temperatures T* only define a S-deconvolution if they are 
an approximate solution of the linear Eqs. (|6|). Only in this case, the temperatures 
T* are similar to the true temperatures averaged by the beam and, consequently, the 
resulting maps contain the right spectra. This fact strongly restrict the methods for 
S-deconvolution. Two of them are described in next section. These equations can be 



written in the matrix form T = AT*, where T and T* are arrays of N x N numbers 
(one number for each pixel) and A is a N 2 x iV 2 matrix. The element A ai weights 
the contribution of the pixel i to the smoothed temperature at pixel a. Quantity A a i 
is assumed to be significant only in the q x q pixels mentioned in § 3.1. 

3.4 METHODS AND RESULTS 

Two methods are used to estimate the S-deconvolved temperature T* corresponding 
to pixel i: In the first one, Eqs. (](]) are solved as a linear system of algebraic equations 
(hereafter LS-deconvolution) where the independent terms are the observed temper- 
atures T a . In the second method, equation (H) is considered as a convolution and, 
then, Fourier transform (FT) and the deconvolution theorem are used to get T* on 
the nodes of the 2D Fourier grid (hereafter FTS-deconvolution). 

3.4.1 LS-DECONVOLUTION 

Iterative methods (Golub & van Loan, 1989; Young, 1971) can be used to solve the 
system (§). In any of these methods, the matrix A is split as follows A=M-P, where 
M is any matrix which can be easily inverted. In matrix form, the iteration scheme 
reads as follows: 

T *(n+1) = M -lp T *{n) + M -l T ^ (g) 

where the superscript n stands for the n-th iteration. The necessary and sufficient 
condition for convergence is that the spectral radius of the matrix Q = M~ 1 P is 
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smaller than unity. This radius is the maximum of | Aj| , where Aj are the eigenvalues 
of Q. Hereafter, the so-called Jacobi method is used. This method corresponds to a 
particular choice of the matrix M. This matrix is assumed to be the matrix formed 
by the diagonal of A, which is denoted A D . A sufficient condition for the convergence 
of the Jacobi method is that matrix A is diagonal dominant (\Aa\ > \Aij\ for 
any i). The dimension of the matrices A, A D and P are iV 2 x N 2 , where N is the 
number of pixels per edge in the map. Since this number is greater than 10 2 in all the 
practical cases, the dimension of the above matrices is very great and they cannot be 
stored. Fortunately, this storage is not necessary; in fact, if Eq. (§) is rewritten using 
indices 

AiiT *(n+i) =Ti _Y^ A l0 Tf n) - £ A l0 Tf n) , (9) 

j<i j>i 

we see that all the elements of the matrix A necessary to get 7i*( n+1 ) can De obtained 
when they are necessary, without storage. Of course, a given element can be calculated 
various times, but no storage is necessary at all. 

Equation (|8]) can be considered as a matrix equation in which each element is a 
block. Matrix A is split in N x N blocks and vectors T and T* in N arrays (blocks) 
of dimension N. Then the resulting N x N blocks appear to have a q-diagonal 
structure and diagonal domination gives limit values for quantities A = e~(A 2 /2cx 2 ) 
and 8 FWHM /A. The maximum value of A and the corresponding maximum value of 
the ratio 6 FWHM /A are given in the first and second columns of Table 1, respectively, 
for various q values. The third column gives the ratio S q /S B , where S q is the area 
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covered by the q x q significant pixels and S B is the area of the beam (a circle with 
radius FWHM ). The top panel of Fig. 1 illustrates the case q — 7. Eleven pixels are 
located inside the beam circle (radius FWHM ). 

Taking into account that, in general, the condition A < \ max is only a sufficient 
condition for the convergence of the Jacobi method (not necessary), we have studied 
numerically many cases in which A > \ max with the hope of getting new convergent 
cases. The result is that the Jacobi method has never converged for A > X rnax . 
This result points out that, in practice, the condition A < \ max is also necessary for 
LS-deconvolution. This information suffices for us. A rigorous mathematical study 
about the necessary character of this condition is not appropriate here. From the 
intuitive point of view, the existence -in practice- of a \ max value is an expected 
result, in fact, in the absence of a \ max , it would be possible to assign deconvolved 
temperatures to billions (an arbitrary number) of small pixels placed inside the beam 
and, furthermore, the right spectrum could be recovered up to the spatial scales of 
these small pixels; this would be a nonsense. 

When LS-deconvolution applies, it is an accurate method. In fact, for 9 FWHM = 
8.8' and A = 4.6875', about 24 iterations suffice to get a very good deconvolved map. 
After these iterations, the relation [££(T/" +1) - T^) 2 } 1 ' 2 / N 2 < 1(T 7 is satisfied 
and, consequently, the method is converging towards a certain map. The question 
is: are the numerical iterations converging to a good S-deconvolved map with the 
right spectrum? In order to answer this question, we proceed as follows: (1) the Fast 
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Fourier Transform (FFT) is used to do a simulation (Si) which does not involve either 
beam or noise, (2) a second map (S 2 ) is obtained -from Si~ by using a certain beam for 
smoothing, (3) the map S2 is deconvolved using the Jacobi method to get a new map 
(S3), (4) the three above steps are repeated forty times (see § 3.2) and, (5) the method 
described in the introduction (see also Saez, Holtmann & Smoot 1996 and Saez & 
Arnau 1997) is used to obtain the modified spectrum in three cases: before smoothing 
(from Si maps), after smoothing (from S 2 maps) and, after deconvolution (from S 3 
maps), these spectra are hereafter referred to as Eu, E 2e , and E31, respectively. If 
deconvolution is a good enough S-deconvolution, spectrum E31 should be comparable 
with E u . Whatever the deconvolution method may be, these five steps allow us to 
analize the resulting deconvolved maps. In all the Figures of this paper which show 
the three above spectra, pointed, dashed, and solid lines correspond to Eu, E 2 £, and 
E 3i , respectively. The top right panel of Fig. 2 shows the resulting spectra for the 
LS-deconvolution under consideration (0 FWHM = 8.8' and A = 4.6875'). We see that 
the dotted line (E u ) is almost indistinguishable from the solid one (E^) for £ < 2000. 
This result qualitatively proves the goodness of the iterative LS-deconvolution. In 
order to compare the spectra Eu and E 3 i quantitatively, the following quantities are 
calculated and presented in Table 2 (entries 7 and 8): The mean, Ml, of the quantities 
1(1 + l)Eg x 10 10 corresponding to the spectrum Eu (col. [3]), the mean, M2, of the 
differences E u — E 3i (col. [4]), the mean MA of \E U — E 3i \ (col. [5]), and the typical 
deviation, E, of the differences of column (4) (col. [6]). The above quantities are 
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computed in appropriate ^-intervals (col. [2]). Entries 7 and 8 show values of \M2\, 
MA and £ much smaller than \M1\, which means that the spectra En and E 3 £ are 
very similar in both ^-intervals (40,1000) and (1000,2000). It is also remarkable that 
\M2\ is much smaller than MA, which means that spectrum E 3 £ oscillates around 
spectrum En giving positive and negative values of En — E u which cancel among 
them. 

Table 2 compares other pairs of spectra displayed in Figs. 2 and 3. Column (1) 
gives the Figure and panel where each pair of spectra are displayed. Tables 3 and 4 
have the same structure as Table 2, but they compare pairs of spectra contained in 
Figs 4 and 5, respectively. The interpretation of the data exhibited in these Tables is 
straightforward. For a given entry, the compared spectra are similar if the quantities 
\M2\, MA and S are much smaller than \M1\. Given two entries with similar \M1\ 
values, the smaller the values of \M2\, MA and E, the greater the similarity between 
the spectra (the better the deconvolution if we are comparing En and E^t spectra). 

3.4.2 FTS-DECONVOLUTION 

The method based on the FT only can be used in the case of small enough coverages 
(almost flat regions) allowing a uniform pixelization. The region covered by the 
observations should be a square with no much more than ~ 20° per edge; thus, 
curvature can be neglected and the covered area can be considered as a square where 
the angles 9 and (f) play the role of cartesian coordinates. Equation (§) can be then 
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seen as a convolution of the function U(9,<f>) = T*{6, <fi) sin(6 l ) with the Gaussian 
beam function W = —^e^ 9 ^ e " >2+ ^~^'^ 2 ^ 2a2 , where coordinates 9' and 6' define the 
observation direction n. Then, the deconvolution theorem ensures that the Fourier 
transform of function U is 

U(k) = T(k)/W{k) , (10) 

— * 

where k is a vector in the 2-dimensional Fourier space. Given a smoothed map and 
a window function, we can find their Fourier transforms T(k) and W(k) and, then, 
Eq. (pT0| ) plus an inverse FT allows us to find the deconvolved map. Unfortunately, 
the use of the FT is not compatible with spherically assymmetric rotating beams. 
If one of these beams measures in such a way that its orientation changes from 
measure to measure, Eq. (|3|) is not a convolution anymore and, consequently, the 
FTS-deconvolution does not apply; hence, the FT can be used either in the case of a 
spherically symmetric beam or in the case of a nonspherical nonrotating beam which 
measures preserving its orientation. 

The maximum value of the ratio 9 FWHM /A compatible with FTS-deconvolution 
has been derived using simulations. In all the 20° x 20° simulations, we have taken 
N = 256 (A = 4.6875'), while quantity 9 FWHM has been varied appropriately. The 
code for FTS-deconvolution has been run in each case. This code follows the five 
steps of the process described above for analyzing deconvolved maps. The left panels 
of Fig. 2 show the spectra En, E 2 e, and E 3i , for different values of 9 FWHM . The top 
left panel, which corresponds to 9 FWHM = 10', shows that the spectra E u and E 3 i 
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are very similar in all the ^-interval (2,2000). For 9 FWHM ~ 11' (middle left panel), 
these spectra are very similar for 2 < t < 1000, while they become a little diferent in 
the interval (1000,2000). For FWHM > 11' , the differences between Eu and E 3 g grow 
rapidly and, for FWHM = 11.5' (bottom left panel), these spectra are very different. 
This is quantitatively confirmed by the numbers presented in entries 1 to 6 of Table 
2, where we see that, in the interval (1000,2000), quantities \M1\, MA and E increase 
as FWHM does. The maximum value of 9 FWHM /A appears to be ~ 2.3. This means 
that, in the case 1.87 < 9 FWHM /A < 2.3, the FTS-deconvolution applies and the LS- 
deconvolution does not. The most dense grid compatible with FTS-deconvolution is 
shown in the middle panel of Fig. 1. In this case, around 16 pixels can be placed 
inside the beam circle. 

An asymmetric nonrotating beam of the form ([|) has been also deconvolved for 
various values of the parameters a and FWHM . The middle right panel of Fig. 2 shows 
the results of the FTS-deconvolution for a = 1.29 and FWHM = 7.75'. These results 
are good in the full ^-interval (2,2000) (see entries 9 and 10 of Table 2). This choice of 
the parameters a and a simulates an asymmetric beam whose effective FWHM along 
the 6*-axis (0-axis) is 91 = FWHM /a ~ 6' (9% = a 9 FWHM — 10'). This beam and 
the most dense pixelization allowing its FTS-deconvolution are shown in the bottom 
panel of Fig. 1. Around nine pixels are located inside the beam ellipse. Other a and a 
values have been also considered to conclude that FTS-deconvolution is possible when 
both af f /A and a e 2 ff /A are smaller than ~ 2.3 (this constraint is equivalent to that 
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obtained in the spherically symmetric case). Finally, we have compared the effect of 
a symmetric beam with FWHM = 10' and that of the asymmetric nonrotating beam 
described above. In order to do this comparison we have obtained forty S2 maps with 
each beam. The spectra obtained from these S2 maps are displayed in the bottom 
right panel of Fig. 2 and quantitatively compared in entries 11 and 12 of Table 2. 
The solid (pointed) line corresponds to the asymmetric (symmetric) beam. Results 
show that the deformations of the original spectrum produced by these beams are 
different; namely, that the E 2 e spectra are distinct (significant asymmetry). 

Vanishing instrumental noise has been assumed so far; nevertheles, partial cover- 
age introduces a kind of sky noise (an uncertainty). In order to estimate this noise for 
a coverage of forty 20° x 20° maps, we compare the SI spectrum extracted from these 
maps with the theoretical spectrum (which would correspond to many realizations of 
the full sky). Both spectra are presented in Fig. 3, where the solid (pointed) line 
corresponds to the SI (theoretical) spectrum. The quantitative comparison of these 
spectra is given in entries 13 and 14 of Table 2. We see that the deviations with re- 
spect to the true spectrum -produced by the partial coverage under consideration- are 
greater than those produced by good S-deconvolutions (compare entries 13 and 14 of 
Table 2 with the pairs of entries 1-2, 7-8, and 9-10. Compare also the corresponding 
panels in the Figures). The deviations decrease as the coverage increases. 
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4 BEAM DECONVOLUTION IN NOISY MAPS 

In § 3, negligible uncorrelated noise has been assumed, thus, for a given beam, a min- 
imum pixel size for LS-deconvolution and another one for FTS-deconvolution have 
been found. These minima define theoretical restrictions for admissible pixelization; 
nevertheless, in the presence of uncorrelated noise, stronger restrictions on pixel sizes 
could appear and, consequently, S-deconvolution could be impossible for some sizes 
close to the minimum size obtained in the absence of noise. Which is the minimum 
size allowing FTS-deconvolution in the presence of a certain level of uncorrelated 
noise? We are going to study this question. 

As it is well known, the amount of noise in a map depends on the observing 
time per pixel, t pix , which is inversely proportional to the pixel area. At pixel i, 
the noise contributes to the temperature an amount 5T* . It is assumed that the 
noise is uncorrelated and has uniform variance a 2 N ] i.e., {5T^ 5Tj) = u^Sij. The 
relation a N = s/(t pix ) 1 / 2 can be used to estimate the level of uncorrelated noise in 
the pixelized map, where s is the detector sensitivity (see Knox, 1995). For s = 
200 nKyf(sec) , A = 4.6875', and a year of uniform full-sky coverage, the level of 
noise is a N ~ 93.9 fiK. Furthermore, using the C t numbers of § 2 with A = 4.6875' 
and 6 FWHM = 8.8', the expected signal S = [(1/4tt)(£/(2^ + l^e"^ 2 )] 1 / 2 takes 
on the value S ~ IOS/jK. Therefore, for the above choice of s, FWHM , and A, the 
signal to noise ratio S/a N is close to 1 and, consequently, noise cannot be neglected 
a priori in order to do beam S-deconvolution. In other realistic cases, the situation 
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is similar. Since a N is proportional to (tpi^ 1 / 2 , too long observation times would 
be necessary to rise S/a N significantly. Fortunately, technological progress leads 
to smaller and smaller s values and, accordingly, the ratio S/a N increases. In the 
PLANCK project of the European Spatial Agency there are two instruments: the Low 
Frequency Instrument (using radiometers) and the High Frequency Instrument (using 
bolometers). One of the radiometers [y = 100GHz and Ofwhm = 10'), working at 
20° K, would produce a noise a N ~ 54 u,K (for pixels with A = 4.6875') during 
a year of uniform coverage; moreover, for one of the bolometers [y = 143GHz and 
Ofwhm = 10.3'), which would work at ~ 0.1° K, the noise would be a N ~ 16.4 \xK 
for the same pixels and time coverage. These data are taken into account below in 
order to analyze the perspectives of beam S-deconvolution in the framework of the 
most accurate project for anisotropy detection in small angular scales (the PLANCK 
mission of the European Spatial Agency). 

In spite of the fact that LS-deconvolution has been very useful in order to analyze 
and understand the existence of a minimum size for pixelization, in practice, only the 
FTS-deconvolution has been used -so far- in the noisy case. The maximum value of 
9 FWHM /A compatible with FTS-deconvolution depends on the level of uncorrelated 
noise. For this type of noise, if the average in Eq. (Q) is performed on many sky 
realizations, the resulting C a (a) values must be very small; nevertheless, if the average 
is done in a 20° x 20° patch, the C a {a) values can be relevant, which means that, 
on the patch, the noise is not properly uncorrelated (its spectrum is unknown). This 
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fact is important in order to understand the effect of this noise on FTS-deconvolution 
of 20° x 20° maps. In the presence of a certain noise which is independent on the 
signal, FTS-deconvolution could be performed using the so-called optimal Wiener 
filter (Press et al. 1988). In such a case, the spectrum of the function U -defined 
above- is estimated as follows: 

U(k) = T(k)$(k)/W(k) , (11) 

where 

m = i Jg@E , . (i2) 

\U(k)\* + \N(k)\ 2 

In the absence of noise, function $ takes on the form $(&) = 1 and Eq. (|TT|) reduces 
to Eq. ([H]). Equations (|TlD and fll2"D cannot be used in practice to deconvolve the 



beam (unknown spectrum of a given 20° x 20° noise realization); nevertheless, these 
equations are useful to understand why the noise can be neglected in some cases. The 
maximum of the \U(k)\ 2 values corresponding to forty 20° x 20° simulations -based on 
the model of § 2- has been estimated to be 1.7, while the maximum of | N( k) | 2 obtained 
from the same number of simulations of pure noise (uncorrelated in great regions) has 
appeared to be proportional to the level of noise a N . For a N = 16.4 fiK, the resulting 
maximum is 2.7 x 1CT 3 . In this case -and also for any current or planned experiment- 
the amplitude corresponding to \N(k)\ 2 is much smaller than that of \U(k)\ 2 . This 
smallness -relative to that of the signal- indicates that, in realistic noisy cases, the 
filter function is close to $(/c) = 1 and the following question arises: Is it possible 
to take $(&) = 1 (noise neglection) to reverse beam smoothing? No theoretical 
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arguments have been found to answer this question. Numerical simulations have 
been necessary. Results obtained from simulations are displayed in Fig. 4. In all the 
cases studied, the noise has been neglected and Eq. (|T0D has been used to perform 
FTS-deconvolution. Good results indicate that noise neglection is appropriate. 

In the left panels of Fig. 4, the level of noise is 16.4 \iK . For 9 FWHM = 8.8' (top left 
panel) spectra E^ and E 3 i are quasi indistinguishable in the ^-interval (2,2000). For 
FWHM ~ 9.5', a small difference between these spectra appears for 1000 < £ < 2000 
(middle left panel). Finally, Eu and E 3 £ are clearly different for FWHM = 10' (bottom 
left panel). Hence, the maximum value of 8 FWHM /A is close to 2 (around 12 pixels 
inside the beam circle). The right panels of Fig. 4 show the same analysis as the 
left panels for a noise level of 54.6 [iK. In this case, spectra Eu and are similar 
in all the interval (2,2000) for 6 FWHM < 8' (top right panel), small discrepancies in 
the interval (1000,2000) have already appeared for FWHM ~ 8.5' (middle right panel) 
and, finally, for 9 FWHM ~ 9' these discrepancies are important (bottom right panel). 
For this level of noise, the maximum value of 8 FWHM /A is close to 1.8 (around 10 
pixels inside the beam circle). This qualitative analysis of Fig. 4 is confirmed by 
Table 3 where quantities Ml, M2, MA, and S are presented in all the cases. 
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5 PIXELIZATION AND DECONVOLUTION IN 
REALISTIC EXPERIMENTS 

Limitations of the S-deconvolution process have been discussed along the paper. If 
the pixel size A is taken to be similar to the beam radius FWHM , the angular power 
spectrum can be only estimated for t < £ max with £ max ~ 180/A ~ 180/9 FWHM ; 
however, for A ~ 9 FWHM /2, the angular power spectrum can be evaluated up to 
(-max ~ 360/^ F1VHM . In these formulae, A and FWHM must be written in degrees. We 
see that, for FWHM ~ 10', the pixelization A ~ 8 FWHM (A ~ 9 FWHM /2) allows us to 
get the spectrum up to £ max ~ 1080 (£ m ax ~ 2160); therefore, for a given beam, the 
choice of the best feasible pixelization is crucial in order to get maximum information 
from observations. The minimum pixel size compatible with S-deconvolution (for the 
methods used in the paper) is hereafter denoted A DE . 

In realistic experiments, various effects -apart from deconvolution- conditionate 
the choice of the most appropriate pixelization. In order to discuss these effects, let us 
focus our attention on PLANCK project (see Tauber, 1999). With a telecope having a 
diameter D, the minimum pixel size allowed by difraction is roughly A DI ~ 1.22c/ Dv, 
where c is the speed of light and v its frequency; hence, this minimum size depends on 
v. In a multifrequency experiment, there is a minimum pixel size A DI corresponding 
to each frequency; for example, in the PLANCK mission (D = 1.5 m), observations 
will be carried out in nine different frequencies ranging from 30 GHz to 857 GHz 
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and, consequently, the size A DI ranges from 30' to 1'. Furthermore, in the PLANCK 
case, the line of sight will move on a big circle in the sky each minute; hence, if two 
successive temperature assignations on the circle are performed at an angular distance 
Aa, the time -in seconds- between these asignations is At = 2.78 x 10~ 3 Aa. The 
angle Aa must be chosen in such a way that (1) no large overlaping of contiguous beam 
positions occurs and (2) time At is greater that the response time of the bolometers. 
For the chosen period of one minute and Aa > 29 FWHM , there is no overlaping and 
condition (2) is satisfied for the PLANCK bolometers. A certain pixel size is only 
admissible if technology plus observational strategy ensure that each pixel is observed 
a large enough number of times during the mission. Let us estimate this number for 
PLANCK. For a pixel size A, the total number of pixels is N = 1.5 x 10 8 A -2 and 
admitting uniform coverage during a year (for qualitative estimates) each pixel is 
observed for a time At p = 0.2A 2 s; therefore, the number of observations per pixel 
is N p = At p / At = 72A 2 / Aa. Finally, for Aa = 29 FWHM , one easily see that the size 
necessary to obtain N p observations by pixel (during a year of PLANCK mission) is 
A* = \{N p 6 FWHM ) l l 2 . Since the FWHM values for PLANCK detectors range from 
~ 30' to ~ 4.5', assuming N p > 100, we see that A* ranges from ~ 9.5' to ~ 1.3'. 
This means that, in order to have a number of observations by pixel greater than 100, 
the pixel size must be greater than A*(N P = 100) = Aioo- 

Given a frequency, there is an optimum pixel size, A OP , which will be assumed 
to be the maximum of the three above sizes Ade, Adi and Aioo- The value of A OP 
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depends on the frequency. 

In order to separate the foregrounds and the cosmological signal in multifrequency 
experiments, various frequencies and a unique pixelization must be used and, con- 
sequently, the best pixelization would be the maximum of the A op optimal sizes 
corresponding to the involved frequencies. This maximum corresponds to the lowest 
frequency under consideration; for example, in the PLANCK case, if all the fre- 
quencies from 30 GHz to 857 GHz are considered, the minimun admissible pixel 
appears to have a size A = A DI (v = 30) ~ 30', for which, the spectrum can be 
only estimated up to £ max = 360. Of course, we could consider only the frequencies 
greater than 53 GHz (with some loss of information) and, then, the minimum pixel 
is A = A DI (v = 53) ~ 17' and l max = 630 and so on. 

Let us reconsider the radiometer working at v ~ 100 GHz with FWHM = 10', 
which was projected to be inside PLANCK satellite. For this radiometer one easily 
find A DE ~ 5', A DI = 8.4' and A 100 ^ 3'; hence, A OP = A DI = 8.4' (£ max ~ 1290), for 
this optimum pixelization, the level of noise is 5.3 x 10~ 6 and the angular spectrum 
can be obtained for £ < 1290 (top panel of Fig. 5 and entries 1 and 2 of Table 4). 
We now consider the bolometer working at v ~ 143 GHz with FWHM = 10.3', which 
was also proposed to measure CMB anisotropy from the PLANCK satellite (phase A 
study). For this detector we easily find A DE ~ 5', A DI ~ 5.87' and Ai o — 5'; hence, 
A OP = A DI = 5.87' (£ max ~ 1840). For the pixelization A = 5.87' the level of noise is 
1.28 x 10~ 5 and we have verified that FTS-deconvolution leads to the right spectrum 
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for t < 1840 (bottom panel of Fig. 5 and entries 3 and 4 of Table 4). 

For the bolometer working at 857 GHz with FWHM = 4.4', we get A DI = 1' and 
Aioo — 1', but the estimate of A DE is problematic as a result of the high level of 
noise of this bolometer. Perhaps, in this case, maximum entropy or wavelets could 
give good results; for example, wavelets could be used to lower the noise before beam 
deconvolution. 

6 DISCUSSION AND CONCLUSIONS 

We expect that beam deconvolution will be important in order to study some aspects 
of the observational maps given by experiments as PLANCK. As an example, let us 
argue that the study of the statistical properties of a given observational map should 
be performed after deconvolution. In fact, various methods can be used to know if the 
maps are Gaussian or they obey other statistics; among them, the estimation of the 
correlation function of pixels where the signal is above a certain threshold (excursion 
sets, Kaiser 1984) and the local analysis of the spots distributed in the map (Bond 
and Efstathiou, 1987). Since the beam smoothes the map, it alters the correlations 
of excursion sets and the structure and distribution of the spots; hence, the above 
methods for analyzing statistics should be applied after a good deconvolution. 

The separation of the cosmic signal and the foregrounds requires a unique appro- 
priate pixelization. In the PLANCK case, we have seen that the optimal size for this 
pixelization, A OP , coincides with the size A DI ~ Q FWHM corresponding to the lowest 
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frequency under consideration; nevertheless, other studies can be imagined (statis- 
tical analysis et cetera) which could be performed on the maps corresponding to a 
given frequency (without previous separation). 

For small enough (but feasible) values of a N plus a certain beam (either spheri- 
cal with a FWHM or asymmteric), our codes allow us to find the most appropriate 
pixel size for deconvolution A = A DE . For the corresponding pixelization, FTS- 
deconvolution leads to a good estimation of the angular power spectrum in the most 
wide ^-interval. The size A DE must be compared to Ai o and A DI to choose the most 
appropriate pixelization A = A op . In the case A OP > A DE , the study about beam 
reversion presented in § 3 to § 5 proves that S-deconvolution can be performed using 
very simple methods. For levels of noise much higher that those of previous sections, 
further study is necessary; maximum entropy, wavelets or other methods should be 
tried out. 

The goodness of a certain pixelization against beam S-deconvolution has appeared 
to be weakly dependent on the particular mathematical method used to reverse the 
beam average. This fact suggests that deconvolution procedures different from those 
of this paper could alter its results. Altough this suggestion should be a motivation for 
studying new methods to get approximate solutions of Eqs. (6) (S-deconvolutions), 
the structure of the system of linear equations to be solved is always the same and, 
consequently, all the mathematical methods could exhibit similar limitations to solve 
it. Indeed, we believe that new deconvolution methods could lead to some modifica- 
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tions of the results of this paper, but not to very different values of A DE . 

Our estimates show that, in the absence of any problem with beam asymmetry, 
current technology could lead to a good estimate of the angular power spectrum (of 
the total signal including foregrounds) in a ^-interval which depends on frequency. In 
the particular case of two instruments on board of PLANCK, we have used optimum 
pixelization and forty 20° x 20° maps to find that the spectra are recovered from 
£ = 200 to i ~ 1300 in the case of a radiometer and from i = 200 to i ~ 1800 for 
a certain bolometer. Since the sky can be divided into ~ 100 of these maps and we 
only need about 40 for a good estimate of the spectrum (for large t values), we can 
select the best forty maps; namely, the maps having minimum contaminations. The 
uncertainty produced by this partial coverage appears to be a little greater than the 
errors produced by the implemented deconvolution procedures (this means that these 
procedures are good enough for us). 

The problem with the deviations of the beam structure with respect to spherical 
symmetry deserves much attention. As discussed above, FTS-deconvolution is com- 
patible with beam asymmetry if the beam orientation is preserved from measurement 
to measurement. If the experiment is designed in such a way that the beam does 
not rotate, our codes for FTS-deconvolution work (see § 3.4.2 and the middle and 
bottom right panels of Fig. 2); however, if the beam rotates, FTS-deconvolution 
does not apply and, moreover, operative methods for making beam S-deconvolution 
are not known; hence, if the assymetry is high enough, beam S-deconvolution is not 
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feasible (so far). In short, excepting the case of negligible deviations with respect to 
spherical symmetry, any effort directed to maintain unaltered the beam orientation 
during observations seems to be of great interest. Unfortunately, in spatial projects 
as Planck, the design of the observational strategy does not preserve this orientation. 

Given an asymmetric rotating beam, it would be interesting to study the whole 
effect of asymmetry plus changing orientation. Even if S-deconvolution is not fea- 
sible, the estimation of this whole effect could be a further direct application of the 
techinques used in this paper. The following method seems to be appropriate: (i) 
average the asymmetric beam -on appropriate shells- to get a new associated one 
with spherical symmetry, (ii) simulate forty SI maps, (iii) smooth the SI maps with 
the assymmetric beam taking into account the orientation change produced by the 
observational strategy; thus, we obtain forty S2 maps, (iv) deconvolve the S2 maps 
with the spherically symmetric beam of reference to get the S3 maps, and (v) esti- 
mate the spectra Eu, E 2 e and, E 3 £. If Eu and Eu are similar enough, the assymetry 
can be neglected, on the contrary, the differences between these two spectra can be 
considered as a measure of the whole effect of asymmetry plus rotation. This study 
should be developed for realistic beams and observational strategies, which is out of 
the scope of this paper. 
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FIGURE CAPTIONS 



FIG. 1.— Top (middle) panel shows the relation between the size of a circular beam 
and that of the pixels for the most dense pixelization compatible with LS (FTS) 
deconvolution. Bottom panel shows the same for an asymmetric beam and FTS- 
deconvolution 

FIG. 2.- Each panel shows quantity £(£+l)Et x 10 10 versus \og£. No noise is present 
and A = 4.6875'. All panels, excepting the bottom right one, contain three lines: 
pointed line gives the spectrum before beam smoothing, dashed line is the spectrum 
after smoothing and, solid line corresponds to the S-deconvolved spectrum. Left: 
top, middle and bottom panels correpond to FTS-deconvolutions with beams having 
°fwhm = 10 '> °fwhm = n '> and fw h m = n - 5 '> respectively. Right: top (middle) 
panel shows the same as the left panels for LS-deconvolution and FWHM = 8.8' (for 
FTS-deconvolution and the asymmetric beam defined in the text). The dotted (solid) 
line of the bottom right panel gives the spectrum after smoothing for a spherically 
symmetric beam with a = 10' (for the asymmetric beam of the text). 

FIG. 3.- The quantities represented are the same as in all the panels of Fig. 2. Solid 
(pointed) line is the modified spectrum extracted from a coverage of forty 20° x 20° 
maps (the theoretical modified spectrum). 
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FIG. 4.- The same as in the left panels of Fig. 2. Left: the level of uncorrelated 
noise per pixel is a N = 16.4 fiK. Top, middle and bottom panels correpond to 

°fwhm = 8 - 8 '? d FWHM = 9 - 5 '> and fwhm = 10 '> respectively. Right: the same 
as in left panels for a N = 54.6 \iK. Top, middle and bottom panels correpond to 
°fwhm = 8 '> °fwhm = 8 - 5 '> and 9 fwhm = 9 '> respectively. Pixel size is A = 4.6875' in 
all cases. 

FIG. 5.- The same as in Fig. 4. Top panel shows the results of FTS-deconvolution 
for a radiometer of the PLANCK mission described in the text. The pixel size is 
A = 8.4'. The bottom panel shows the same for a bolometer of PLANCK and 
A = 5.87'. 
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TABLE 1 

CONVERGENCE OF THE 
JACOBI METHOD 



Q ^max {@fwhm I '^)max S q /S g 

3 0.5000 1.998 0.72 

5 0.4565 1.879 2.25 

7 0.4559 1.877 4.43 

9. 0.4559 1.877 7.32 
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COMPARING SPECTRA 
FIGS. 2 AND 3 
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TABLE 3 

COMPARING SPECTRA 
FIG. 4 
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TABLE 4 

COMPARING SPECTRA 
FIG. 5 



Panel ^-Interval Ml Ml MA E 

Top 40-700 8.96 0.12 0.19 0.22 

Top 700-1290 -14.99 -0.22 0.27 0.33 

Bottom 40-1000 9.43 0.13 0.23 0.30 

Bottom 1000-1840 -16.24 0.31 0.40 0.49 
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